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ABSTRACT 



There are many factors, such as aircraft conf iguration 
and weight, winds aloft, airspeeds flown, altitude, distance, 
etc. , which affect fuel constimption in turbojet aircraft o 
for any given combination of these factors a flight path can 
be determined that will result in the least fuel constimed for 
a ground distance covered. Under divert conditions from air^ 
craft carriers at sea to fields ashore the choice of the 
optimal flight path is critical . The many possible combi- 
nations of factors lead to the adoption of computer flight 
planning. Pilots can avail themselves of computer solutions 
during flight planning and briefing sessions, and after 
take-off can receive further information via UHF radio. 
Typical flight handbooks display fuel flow data, etc. in 
such a manner that the pilot must "guesstimate" entry 
parameters such as average horizontal weight, or weight prior 
to descent. Several iterative procedures are developed that 
provide exact solutions to these Important figures. Thus 
the computer flight planning system will provide more 
accurate solutions, and free the pilot from this chore so 
that he may better spend his time briefing tactics. 
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CHAPTER I 



PURPOSE OF THE STUDY 

Naval Aviation has been described as hours and hours 
of boredom interrupted by moments of sheer terror^ The 
simple word •' BINGO " provides for some of these anxious 
moments. Bingo is well known to everyone as a friendly 
game of chance, but the word carries a special meaning to a 
U,S, Navy carrier pilot. It causes concern both ashore and 
at sea, and especially in the cockpit of at least one air- 
plane that is flying circles around some aircraft carrier 
somewhere at sea, 

"portress 501, this is Atlas tower. Bingo 100°/125 
miles. Change to departure control, 316, 6, for radar 
vector, " 

"Atlas tower, this is Fortress 501, Roger, switching 
316,6, " 

This conversation might have occurred west of San 
Clemente Island off the coast of Southern California. To 
the pilot of Fortress 501 the message is unmistakable. He 
is to turn to 100° magnetic and land at the field 125 miles 
away, rather than attempt to land aboard Atlas, U.S.S, 
INTREPID, CVA-11, 

In most cases an aircraft will receive a BINGO because 
of a flight deck accident, or because the pilot is having 
difficulty landing aboard the carrier due to bad weather 
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and/or a badly pitching flight deck. The high fuel con- 
sumption of jet aircraft at low altitudes precludes having 
the pilot hold until the deck is clear, or in the latter 
case attempting more landings aboardo 

liJhen the pilot of Fortress 501 switches to 316® 6, 
he will receive a radar vector to the nearest suitable 
landing field, the weather at the field, and the approach 
control and tower frequencies. But he will have to deter- 
mine his own flight profile; that is, how high to climb, 
what airspeeds to fly, and at what distance from his 
destination to commence an idle descent. Since the air- 
craft is usually low on fuel when the pilot commences the 
BINGO, the pilot attempts to fly the flight profile that 
minimizes fuel burned. There are no other constraints on 
the pilot, for he is only interested in running out of 
ocean before he runs out of fuel. 

The factors that affect fuel consiunption in a jet 
aircraft are altitude, aircraft weight, and of course, the 
air distance travelled. Jet aircraft burn enormous amounts 
of fuel at lovi altitudes and during the climb to higher 
altitudes. During an idle descent from altitude the fuel 
flow is only one-sixth that of the climb portion. A 
heavier aircraft burns more than a lighter one, and a head- 
wind increases air distance and thus Increases fuel burned. 
Jith these factors in mind, the pilot must decide what com- 
bination of climb, cruise, and descent flight paths will 
take him to his destination with the least fuel consumed. 
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A wrong choice could result in the loss of an aircraft and 
perhaps a pilot. A lesser error might result in a safe 
landing, but several hundred pouinds of jet fuel will have 
been wastedo This latter error might also result in some 
premature gray hairs on an intrepid naval aviator. 

Since a determination had to be made at what fuel 
weight to BINGO the pilot, the knowledge that all pilots 
would fly an optimal (minimum fuel) flight path would allow 
the BINGO to be delayed for perhaps one more pass at the 
flight decks Also, any "gravy" could be reduced from the 
bingo weight. This lower weight would increase operational 
readiness by keeping all aircraft aboard the carrier where 
they belong, and where all commanders want them. If an 
aircraft is sent to the beach, the carrier must remain in 
the area to send messages concerning ship's position and 
overhead times so that the stray bird can come home to 
roost. Since carrier skippers like to hide their ships 
out at sea, there are obvious tactical advantages in 
keeping all aircraft aboard. 

An ideal way to ensure that all pilots fly an optimal 
flight path for all distances, aircraft weight, and wind 
conditions would be by computer flight planning. Thus when 
Portress 50 I contacted departure control on }l6 ^ 6 he might 
hear, 

"Roger, Portress 5 OI 0 Climb at 285 knots to 19*000 
feet. Cruise at 260 knots indicated. Commence idle descent 
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when 44 miles from destination flXo" Departure control could 
then follow with the standard weather and radio aids infor^ 
mationo 

This thesis considers only the short BINGO problem in 
its development of a computer flight planning program. This 
applies to a small, but important segment of aviation. The 
obvious extension is to use this type of flight planning for 
all long distance cross-country flights. The only modifi- 
cation necessary is to incorporate several horizontal legs 
at altitude sandwiched between a climb and a descent leg, 
rather than just one as does the BINGO program. The cross- 
country program would also require a change of optimal 
flight altitude to a higher one when the aircraft weight 
decreases enough to warrant a climb in order to remain 
optimal. 
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CHAPTER II 



REPORT OP THE STUDY 



The BINGO problem may be stated mathematically as a 
classic minimization problem: 



min 
S « t e 



(“'o-'S* 



h<40 X 10''^, in even thousands 






> B 



Where B = aircraft empty weight 
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Pigure 1 

Definitions and abbreviations: 

Nq = aircraft weight at start of flight 
= weight at end of climb portion 
/l^ - weight at end of horizontal portion 
= weight at destination 
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h = altitude in thousands of feet 

X = ground covered in climb portion 
c 

= ground covered in horizontal portion 

X = groimd covered during idle descent 
d 

X = total distance to destination 
t 

t = time 

CAS = calibrated airspeed, velocity of the air over 
the pitot tube 

TAS = true airspeed, CAS corrected to sea level 



TAS = 




where 




= density of the air at altitude h 

Pq = density of air at sea level 

GSPD = ground speed, aircraft speed relative to the 
ground (TAS corrected for wind) 

F = fuel burned in pounds 

1/f = rate of fuel flow In pounds per nautical mile 

DI = drag index for the particular aircraft con- 
figuration 

A flight profile depicted by Figure (1) is used because 
of standard operating procedure, and also because the air- 
craft Instrumentation is such that it is the only profile 
that a pilot can fly accurately An arcing path such as a 
semi-circle or a cycloid is impossible* 

In order to fly a flight profile in the form of a 
cycloid or a semi-circle a pilot would have to rely on his 
VSI (vertical speed indicator) to accurately control his 
rate of climb throughout the entire flight. This instrument 
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is inaccurate at best since it is a pressure instrument and 
fluctuates with slight changes in pressure* The VSI also 
tends to lag the true pressure changes* A climb* level* 
descent flight profile can be accurately flown because the 
pilot need only control his airspeed and altitude* Both 
these cockpit indications are excellent in all aircraft, and 
a pilot gets acquainted with them from his first step into a 
cockpit. 

The objective function* is the fuel consumed 

in travelling the distance X^, since the only weight loss 
will be the fuel burned. The constraint > B is obvious. 

If the destination cannot be reached, even by flying an 
optimal flight profile, then the aircraft should rendezvous 
with a tanker aircraft if one is available, or else continue 
landing attempts aboard the carrier. It is not obvious that 
X^ + < X^ - X^ is a necessary constraint, A possible 

flight path could be to continue horizontally until the 
destination is reached, and then commence a circling descent 
such that X^ = 0, Navy flight tests show that an idle 
descent at the proper airspeed for X^^ will use less fuel 
than any other possible flight path covering the distance 
X^« This type of idle descent is standard operating proce- 
dure for jet aircraft. 

The altitude, h, must necessarily be less than or 
equal to the service celling of the aircraft being flown. 

For this study h < 40,000 ft,, the service celling for the 
A4c, Is used. Even thousands of feet are used, resulting 
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in ^0 possible flight profiles o Using every 500 feet will 
result in 80 possible profiles, and using every 250 feet 
will result in l60„ This study uses even thousands merely 
for ease of presentation^ Only a minor change is required 
to solve for any desired niimber of possible profiles® 

The following assumptions were made concerning the 
flight path of Figure 1® 

1® the aircraft makes an Instantaneous transition 
from the climb to the level flight attitude, and 
from the level to the descent® 

2® all airspeed changes are instantaneous® 

3® the aircraft can remain on track despite any 
cross-wind® 

4® wind information is known, and there are no up- 
drafts or downdrafts® 

performance data such as fuel flow and airspeed were 
obtained from NATOPS Flight Manual for the a 4C aircraft. 
Although only the a4c Is considered here, similar data can 
be used for other aircraft in the inventory® All perform- 
ance data was gathered by U®S® Navy flight tests® The 
peculiar form of the graphs (Figure 6=13) makes it diffi- 
cult to form mathematical functions for such values as 
maximum range airspeed,, This precludes formulating the 
problem as a standard Lagrange minimization problem® 

The computer program developed is one that calculates 
the fuel burned for all altitudes up to 40,000 feet, the 
service ceiling of the a4c, and chooses that flight profile 
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that results in the least fuel burned,, The first altitude 



tried Is level at 1000 feet for the entire distance, X o 

t 



The remaining thirty-nine are flown as follows c 

(1) a 100^ power climb to h- thousand feet 

(2) level at h-thousand feet at the max-range air- 
speed, (that airspeed that gives the most miles 
per pound of fuel ) 

(3) an idle descent at that airspeed that covers 

with the least fuel burned 
d 

Aircraft performance data is presented as follows® 

CLIMB PORTION; 

P = P (Wq, h, DI) 

CAS = C (i , h, DI) 

u 

time = t h, DI) 

HORIZONTAL PORTION: 



f 



P 




h, DI) 



( 1 ) 



P 



P(X^j -g-- ' , h, DI ) 



(2) 



CAS = . h, DI) 



DESCENT PORTION 



P = PCNg, h, DI) 

CAS = CCN^s DI) 

Xd = XCWg, h, DI) 

time = t(h, DI) 



( 4 ) 

( 5 ) 



( 3 ) 
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Future functional notation will omit the drag index, Dio 
It will be held constant at DI - 50, the value for an air- 
craft with no external stores o This is almost always the 
case in a BINGO situationo 

The first flight profile tried is level at h = 1000 
feet. Since this differs from the others only insofar as 
there are no climb or descent portions, a typical climb, 
level, descent profile will adequately explain both possi- 
bilitieso For the climb to any altitude, Wq and h are knowno 
Prom Figures 6, 7, and 8, F, CAS, and the time to climb 
are easily determined. CAS is converted to TAS, and the 
predicted wind information (averaged over h-thousand feet) 

converts TAS to GSPDo - F = W. , and GSPD x time = X o 

0 1’ c 

Each 5000 feet the CAS is changed to the optimal climb 

speed for the next 5000 foot portlono If h > 5 the climb 

portion of the flight results in ([[h/5] + 1) iterations. 

At the start of the horizontal portion, h, and 

(X “ X ) are knoTm. Figures 9 and 10 however, require the 
U Q 

average weight during the horizontal leg, (V/^ + W^)/2, as 
an entry parameter. Recall that the fuel burned during 
the horizontal leg is an implicit function of and W^. 

From equation 2, 

N-, + No 

F = F(X^, , h) 

h 2 

The procedure for using the graph in Figure 10 is to 
enter with the average weight, proceed horizontally to the 
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altitude h, read vertically downward to the DI = 50 line 
and horizontally to the left to read f, the fuel flow® 

Three different methods were tried to solve the im 
plicit function for the average weighty 

METHOD I 

Approximate the altitude lines in Figure 10 vjith 
straight lines as shown in Figure 2. 




Figure 2 
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Prom geometry we obtains 



tan 0 = 






( ^ 3 ^ - 10 , 



t 

000) W-j "I- Wp 

(Ml - 



t = 






) . tan 0 



Knowing t, we solve for fuel flow, f, using 



f = f + t » tan a 
o 

W-t + Wp 

f = f (W - — = -) tan 0 tan a 

o i 2 



Since p - = j} = fi from equation 1 

■p ± ^ 







( 6 ) 



((/here is defined in Figure 2< 



Solving this quadratic for will allow a good approximation 

i/'/t vio 

to the average weight, This approximation is 

believed to be quite accurate, since the altitude lines are 
fairly straight. An exact approach is given In method two. 



METHOD II : 

Allowing for the fact that the altitude lines are non- 
linear, a more precise but lengthy method for solving for 
the average weight is the method of successive approximations. 
Referring again to the graph of Figure 10, we see that the 
average weight is a number such that, if it is used to enter 
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the table, it produces a fuel used, F, such that - F= 



j^hen lAfg is combined with to get the average ( 



Wl ^2 



), this 



average will exactly equal the average we used to enter the 
table. Because of the construction of the table, the true 
average is the only entry number that will give the average 
back again. 

The first approximation to the average is i . The 

1 

aircraft continuously loses weight, thus we know that this 
approximation is too high. Entering the table with yields 



“ F-, ) , 

^ ± ±— ^ . Decrease 



an F^» This yields an average ^ , 

by some l>0o Entering the table with - A. gives F^« 
Form another average: 



= ^ ^ 

2 1 



Continue in this manner until 






F ) 



W- 



=• nX 



This method is easily suited for computations on a digital 
computer, but it proved too lengthy. The greater the accur- 
acy desired, the smaller K must be. But a small K requires 
a large number of iterations. This method was abandoned 
for the third and final method. 



METHOD III: 

Recall from Figure 10, F= F(-i — — , h)» The graph 

2 

shows that the fuel burned is directly proportional to the 
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aircraft weight As a first approximation to the true 
average weight, use „ This is too largOj therefore it 
results in a fuel burned, F^, that Is too large. 




F^ = F(AV^, h) 

From equation 6 

p „ \ 

1 f + (iif - AV- ) tan 6 tan a 
oil 

Let 

= D a constant 

f + tan e tan a = B a constant 
o i 

tan e tan a = C a constant 



The first approximation to the fuel used on the horizontal 
leg is; 



p _ where X^ = AV, 

1 B + CX^ 1 1 



(7) 



Since F^^ is too large, as a next approximation to the average 

weight decrease AV, by g-P, o Thus in equation 7, multiply by 

1 1 

” g, and add to both sides, ie get 
1 



AV = If 
2 1^1 



D 

2(B+ '“^1 
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AV = X = _L_ -r 
2 2 B -i- CX^ 1 



(where A = - 



Entering Figure 10 with X^ yields 



D 



Fp = 

2 B + CX, 



and 



AV_ = X_ = + 

3 3 B + CX^ 1 



This Iterative process yields a sequence 



X = 



X„ = 



A 



B+ CX-, 






Xo = 



• + 

B + CX„ 1 



X = 

B + CX^ , 1 

n - 1 



where 



A 

B = 
C = 



-150 

»34 

-6^4 X 10 
4 

1.5 X 10^ 
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Solvinti for tns first few points yields a nest of closed 
intervals as shown in Figure 3<» 



^2 ^4 



_1 1 1 1 _ 

X6 x^ 



Figure 3 

At this point it can be seen v;hy - §P ) v/as used 

1 

as the second approximation to the average x\reighto Since 
we know that is well above the average, it is desired to| 
decrease it by a quantity large enough to drive the second 
approximation below the average* A number like *8 or *9 
could have been used to multiply F^, but the use of ^ is 
sufficient and causes the sequence to converge faster than 
numbers like *8 or .9. 

The manner in which the successive points alternate 
led to applying the theory of continued fractions in an 
attempt to show that the sequence { kh ^ converges. Hall and 

r2] ^ 

Knight*- show that each successive convergent of a con- 
tinued fraction is alternately less than and greater than 
the true value of the continued fractlono 
Rewriting equation 8 yields 



X = W + __£L 
1 B + CX 



n 
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Back-substituting for ^ by X^ _ ^ - 2 ^n - 3’ 

etc. yields: 



Xj ^ ‘‘^1 



A 

B + C ( + 



A 



B + C ( + 



B *(• C ( "h 



e 



+ 



A 

B+ CX 

o 



The form of the classic continued fraction is obvious. 
Factoring out the A in the numerator, and multiplying in 
each convergent by C yields: 




+ A 






1 



( B + C ) “i- 



CA 



(B+ ) + 



CA 

( B + C ) '¥ 



Q 



Dividing numerator and denominator of each successive 
fraction by (CA) 



r* 



X ■•= + A 



(B+ CIJ ) -r— 

1 /B't-Cv/^ 

”CA~ 



(B + C:i ) + 

1 /B+Cbh 



CA 



In the continued fraction in the brackets, let 
a 

B '“I" 



B -i- Ci'I >0 
1 



CA 



> 0 
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Then 



. * 



e = 



(9 



a + 



t) ''i^ 



a 



b -ir 



As we let replace the 4^^, etCo convergents by 9 



. ^ 



1 



a 



= 



b “i- 0 



a( b 4 0 ) + 1 
b + 0^ 



b + 0 



aih-irB^) +1 



Solving the resulting quadratic yields s 



* - (a^b^4 4ab)^ 

2a 



0* = 



Thus we have that 



X* » u, + A I 

1 \ 2a 



The minus sign is selected because < A 

1 

This third method gives the average weight for any 
given value of and ho For most values ^ the sequenc 

converges In 6 uo 8 iterations, whereas the second method 
required ivore than 50. 
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It is interesting to note that in the continued frac- 
tion, equation 9» for real but unequal odd and even con- 
vergents, Van Vleck's theorem can be appliedo For a or b 
imaginary, Stieltjes theorem applies® Proofs of these 

[ 3 ] 

theorems are given in ,'^all‘“ ® Having solved the implicit 

function for ^ > Figure 9 yields GAS® It is converted 

to TAS, and the X'rind at altitude h changes TAS to GAPD 



time = / GSPD 

air distance = X 



a 



time X TAS 



v^i + ^2 



Figure 10 gives F as a function of X and 

3 . 

Since optimal CAS is a function of the aircraft weight, 
as the weight decreases during the horizontal leg the GAS 
should be changed to remain optimal. In this program, when 
the aircraft weight decreases by 500 pounds, CAS is recom- 
puted. This results in breaking up X^ into segments of 
length Xj^ = 500 X f® This results in ([Xj^/X^] + 1) iter- 
ations of the horizontal leg computations® 

At this point it should be noted that in order to 
solve for the fuel burned in the horizontal leg by either 
of the three methods, we require Xj^. At the start of the 

leg we know only X,. and X . Since X, = X. - (X. -t-X ) the 

U C xl u c cl 

value of X^ must be determined prior to the start of 

the horizontal leg® From equation 4 recall that 

X, = X (>/ , h). The problem resolves to this: in order 

u 2 

to solve for vJ2, we need X,® But we need before we can 

d 2 
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solve for X . Another method of successive approx Imac ions 
d 

is developed that results in a recurrent sequence that is 

shovm to converge 6 Once It Is shown that v;e can find 

(and consequently Xj^,J at the start of the horizontal legj 

then all previous procedures are justifledo 

The air distance covered during the descent Is 

X = X (he, DI, wind)o Prom equation 3 ^ CAS = C Dl)o 

ad ^ 

Prom equation 5* t = t (h^ DI)o 
Prom Pigure 11 




Pigure 4 

CAS ^ .190 + tan (^2 12000) 

Let (I90 ^ 12000 tan 3) = a 
tan 3 - t> 

CAS — a ’t" 

2 



26 



and 



TAS = 9 CAS = k c CAS vjhere k = -i- 

d o o d 

X 9 X, 

time in descent = t , .= - — — - 

^ TAS GSPD 

= t^ X GSPD 
d d 

GSPD = TAS ± wind = k^ x TAS 



where k^ is a constant to correct TAS for wind„ 



GSPD = k- (k X CAS) 
j- o 

G3PD = k (k [a + b^J 1 
1 o'- 2-^ 



Since 



X, = t X GSPD 
d d 

X = t (ak k + bk-k W ) 
d d ' 1 o 1 o 2^ 



Let 



ak^ k^t , = a 

1 o d 

bk^ k t =3 

1. o d 



X = tt 3 If/ 
d 2 



( 10 ) 



Begin the Iterative process by letting (i/ ) = 

2 0 

large value will yield an X, that is too large, 

d 



vi 



This 



Since 



X, - X - (X + X , ) , the resulting X, will be too smallo 
n t c d h 

Fuel burned on the horizontal leg is directly proportional 



to the distance, therefore F will be smaller than the true 

h 
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value., (iV ) = i - P „ This value is still too large, bujf 

2 1 1 h 



is closer to the true value of than was the first approx 



mat ion, 



From equation 6 we have 



X. 



J, - i = 



n 



^ ^ f + VJ. tan 9 tan a - tan 9 tan a ® / *^1 '*■ 

o 1 —T“ 



Let 



f + tan 9 tan a = b 
o 1 



and 



tan 9 tan a 



i. - VJ, 



Xh 



i 2 B+C + 






Solving the quadratic for yields 



-B± (B^ - 4C + X. J ) 



2C 



But from equation (10) 



X - a + |3(i<V ) 
d 2 o 



Xj^ (X^ ^ X^ ) » X^ - (X^ = X_ ) - a - 3(W 



t c 



t c 



2 o 



Let 



(X^ “X - a 



X ^ d ^ 
a o 
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and 



(Wz)^ = 



= B ± (B^ - 4c [-CV/^ BW^ + [d - y (W 2 )q]] ) 
»_ _ 2C 



Grouping constants , 



(^2)1 • A - (B-C(.J 3 )p 



where 



A = 5o 65 X 10^ 
B = l»8l X 10^ 
C = 2.35 X 10^ 



(Wp) 

^ o 



= lo 6 X 10 



The iterative process thus described results in the sequence 



m ■■ 



X 






X, = A" (B- CX 
J- o ^ 

X^ - A - (B- CX^)^ 



Y = A - (B» CX ) 
^ n - 1 



(lOol) 



Calculating the first fevj values results in a monotonic non- 
increasing sequence of points as in Figure 5® In order to 
show that such a sequence converges* we refer to the theory 
of contraction mappings as illustrated by Lyusternik and 

f4] 

Yanpolskii" . 



29 



x^ x^ x^ 



X. 



'0 



Figure 5 

Let f be a continuous operator from £ into £ , and 

n n 

Y = f (X)o In order to solve for the equation X = f (X) 

vre set up the iterative sequence of elements: x^,x^ jXg , « . . » 

where x^ is an arbitrary element of £ , Here we have 
0 n 

= f (X )« If the sequence { Xj^\ is convergent to some 
X'% then X‘"' is a solution of the equation X=f (X), and X^ 
is called a fixed point of the transformation Y=f (X)o 

The principle of contraction mappings provides a con- 
dition for the existence of a fixed point of the transfor- 
mation Y = f (X). 

D£FINITI0NS : 

(1) An operator f from £ into E is said to be a 

n n 

contraction mapping if there exists a constant q,(0<q<l), 
such that for any x* x £ E^ 



f(x ) - f(x)| < q I 



n 



x^-x 



(2) A sequence ^X^^ having the property that, given 

any £>0, there exists a subscript N such that I X - X I <F fo 

' m n' 

all m>N, n>N, is said to be fundamentals 
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BOLZANO- CAUCHY CRITERION; 



A necessary and sufficient condition for the sequence 
^^n^ to be convergent is that it be fundamentals 

The folloi'/lng shows that the principle of contraction 
mappings provides for a fixed point of the transformation 

y = f(X)o 

Let f be a continuous operator from E into E^, such 

2 

that f is a contraction mapping and in C , ioO. twice con- 
tinuously differentiable in all components c Then there 
exists a solution X* to the equation X = f(X)» and this 
solution is uniqueo The iterative sequence formed by 
successive approximations is convergent to X* whatever the 
initial approximation Xq® 



PROOF : 

Form the iterative sequence of elements 



X 

o 



X 



2 



= 

= f(X^) 
= f(xp 



X 



n 



f (X T ) 
n - 1 






( 11 ) 

( 12 ) 



Subtract equation (11) from (12) and take absolute values 



But f is a contraction mapping by hypothesis* 



+ 1 



I f (X ) - f (X ) I < q I X - X I 
n n-ii— In n-1' 



n+1 nl- In n-1' 



( 13 )' 



Xp - X I < q I X - X 
^ 1 ' “ I 1 c 



X - X I < q Xp - X I = q"^ I X - X 
J , ^ ^ 1 1 o 



X ^ - X ' < q 
n + 1 n — 



n 



X - X 
1 o 



(14) 



If the sequence, equation (14), can be shown to be funda- 
mental, then by the Bolzano -Cauchy criterion it is converge: 
Choose ^>0, and N large* To show that 1 X„ - X ' < F for al 

I in ^ >- 



m > N, n>N, we add and subtract all values of X, , m < k < n. 






n+1 n 



1 


- X 


•h 1 X 


- X 1 




1 m 


m-1 


1 m ~], 


m-2 1 


... + j 



From equation (14), 



X - X 
m m-1 



^ m-1 , , 

< q I X. - X 



X , - X ^ 
m-1 ni=_2 



n+1 n 



m-2 I 

< q Mi " ^ 

I 1 o 



n 



< q X - X 
I 1 o 
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j Xj^~X T T ~^,«l o o -K5_^ ^ )*|X-| “X 



m-1 



n-1 n-1 n' - 



o 



< q^(l+q+q^+- ° )»| X^ - X | 



Since q < 1, the partial geometric series is less than 



1 -q 



X - X +X -ooo -X + X ^ X 
Di iu-1 m-1 n-1 n-1 n 



<q" (^) . 

1 - q 



X -X 
1 o 



iAfe can choose n large enough so that the sequence is funda- 
mental and therefore convergent,, Since q<l, we can make 
the right hand side arbitrarily small* In fact, | X^^ - X^ ^j- 
as n-»”«>^ 

From equation (11) 



X = f(X , 
n n-1 



X " f (X ) 
n n 



= I f (X . )-f (X ) 1 <q 'x„ - X I 
' n-1 n ~ I n n-ll 



I X - f (X ) I ^q I X - X I 
In n I - ! n n-1' 



0 



Since f is in C^, as X^-t- X*, f (X^^ )~^f (X*) 



Thus 



and 



X""- f (X"") 



X*^ = f (x'') 



= 0 



Fe D, 
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If I f * (X) I < q < 1, then f is a contraction mapplngo From 
the mean value theorem of differential calculus 

f(b)-f(a) = (b»3.) « f * (X) for some X, a<X<b. 

ft ^ f (b) ° f (a) 
b “ a 

If I f' (x) I < q < 1, then I f (b) = f (a) i < jb - a 1 
and I f(b)=f(a) I < q jb = a | 



From equation 10.1 the 



(f ) = A - (B - CX)‘ 



f* = 



~r 



2(B - CX)2 



function under study is: 



To show that the function is a contraction mapping it is 
necessary to show that jf*| < q < 1 for all values of X 
that vjill be encountered. It is easily seen that when 

r3 

X = ” the derivative is not defined, and that when X > 

C 

It is Imaglnaryo 



The values of the constants in the equation are: 

A = 5c 65 X 10^ 

B = I 08 I X 10^ 

G 2.35 X 10^ 

X = y = 1.6 X 10^ 

O 1 
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o Ito 



z/e desire to show that X < - = 7,7 x 10^ for all X that 

C 

are encountered in the iterative process. The first 
approximation to X is; 

X = = 1,6 X 10^ 

w X 

X^_ = 5o 65 X 10^ - (1,81 X 10^ - 2,35 x 10^x1,6x10^)^ 

X]_ = 5o65 X 10^ - (18,3 X 10®)^ 

X^ = 1,39 X 10^ < X^ 

Xg = 5 o65 X 10^ ^ (1,81 X 10^ - 2,35 X 102 X X^_ 

x^ = Ie35 X 10^ < x^ < Xq 

It can easily be shown that each value of X will be smaller 
than the one that precedes it, since at each iteration a 
larger value will be subtracted from the constant A, Thus 
the sequence X^^ is monotonic non-increasing, and its upper 

D 

bound is vi/ . Since , all values of X that will be 

X X 0 

encountered are also strictly less than — , Solving for 

C 

the first derivative when X = yields: 

I f« I = X 10^ < ^ 

2(4,243 X 10^) 

Thus I f’ I < 1 for all values of X in the iterative sequence. 
This is so because is monotonlc non-increasing, and 

any value of X < increases the denominator of the deriv- 
ative, thereby decreasing its value. Thus f is a contrac- 
tion mapping, and the sequence ^X^"^ converges, Lyusternik 

r4i 

and Yanpol'skll^ shox\r that if f is a contraction mapping, 
then the sequence is convergent to X* as fast as a geometric 
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progression with ratio q. Thus the computer iteration 
procedure is not lengthy since in this case q = « 0027 ® 

The knowledge that we can solve for and thus X^, 
at the start of the horizontal leg of the flight assures u£ 
that we can correctly find the average weight on the leg 
using the procedures described earlier. We can solve 

- F = 1^2' Knowing ^2 we solve for the fuel used in the 
descent using p = f and = ^2 - F, The fuel 

burned during this profile (W^ - W^) is stored, along with 
the altitude and proper airspeeds for the profile. 

Once fuel figures for all profiles are calculated, a 
searching procedure selects the global minimum and prints 
it as the optimal fuel along with corresponding altitude 
and airspeed figures. The possibility of a tie is virtuall 
eliminated by reading fuel figures to four decimal places. 
Should several local minima arise, the search is designed 
so that the global minimum is always selected. 

A listing of the Fortran IV program and a sample out 
put is given in Appendix B. All variable names are written 
such that they are easily recognizable. TAS is true alrspe 
DIST is distance, etc. As inputs, the program only require 
the course and distance to the field, the type aircraft 
(a 4 =1, A? = 2, etc. This allows for the difference in 
aircraft), the weight of fuel aboard the aircraft at time o 
BINGO, and the wind information. The wind can be up-dated 
at intervals as set by meterologica.1 readings. 
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CHAPTfiR III 



RESULTS AND CONCLUSIONS 

For a fixed distance, nautical miles, (100 < < 250), 

the winds aloft and the aircraft weight combine to determine 

the optimal flight profile. Under a zero wind condition 

and normal BINGO fuel weight, the optimal flight profile is 

to climb to the service ceiling and commence an idle descent 

at X,, or to climb until X, is reached, and then to start 
d d 

the idle descent such that X. is zero. For fuel weights 

n 

higher than the normal BINGO weights, the optimal altitude 
is lower than the service ceiling because of the excess 
fuel used in climbing to high altitudes when the aircraft 
is heavy. 

Since zero wind conditions are rarely if ever en- 
countered, the wind velocity aloft is a most Important factor 
in determining the optimal flight profile. In order to 
test the effect of winds on a standard BINGO problem, a 
typical wind pattern for winter months along the coastal 
region of Southern California was obtained from the weather 
facility at the Naval Auxiliary Landing Field in Monterey. 

The following winds are typical of a northern hemis- 
phere cyclonic low pressure area. 
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ALTITUDE 



/fIND 



0 - 5000 ft 
6 -10000 ft 
11 -15000 ft 
16 -20000 ft 
21 -25000 ft 
26 -30000 ft 
31 -35000 ft 
36 -40000 ft 



270 *^/ 10 kts 
240 °/ 15 kts 
210 °/ 20 kts 
180 °/ 30 kts 
150°/ 40 kts 
1200/ 60 kts 
095°/ 80 kts 
0800/ 100 kts 



The BINGO situation considered is such that a pilot leaves 
a carrier off the coast of Southern California and flies 
090° / 160 miles to NAS North Island near San Diego, The 

I 

fuel weight at start of BINGO is 2200 #. | 

For this example a flight profile with the service 
celling as the optimum altitude would require 1499 . A 
guess by a "seat of the pants" acquaintance of mine would 
use 25000 ft as the optimum altitude and would require 
1^35.3#. The computer flight plan predicts the optimal 
flight profile as follows. 

Optimum altitude is 20,000 ft 
Climb speed is 3 l 4, 4 kts 

Cruise speed is 25 o . 0 kts 

Start descent when ' 32 . 1 ' miles out 

Descent speed is 195 . o kts 

Fuel required is 1135 . 6 lbs 

The savings in jet fuel of 300 # for the optimal over i 
the guess, and 364 # for the optimal over the climb to servi 
ceiling is substantial when it is considered that this 
savings would allow two more landing attempts either at the 
carrier or at the field ashore. Multiply this savings by 
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thousands of flights per month if the system were Incorpor-* 
ated in the fleet, and it is easily seen that the jet fuel 
saved would soon become signiflcanto The two extra landing 
attempts would probably result in fewer aircraft ramaining 
on the beach over nighty 

In order to demonstrate the usefulness of extending 
the present BINGO computer program to long distance flights, 
a flight of 600 miles was flovm with a fuel weight of 6 OOO 
pounds at the start 0 The entire distance is flown with a 
course of 115° magnetiCo This is not realistic since even 
Nest - East cross-country flights sometimes require course 
changes of more than 20 ^ « Also, the present program does 
not change altitudes to higher ones when the aircraft weight 
decreases to a value that suggests a climb in order to 
remain optimal « Nevertheless, some fuel values for 
different profiles indicate the savings in fuel that would 
be realized if the optimal profile is floxvno For example: 



Fuel (Ibs^ ) 



Altitude (ft. ) 




10,000 



5>000 



4333.9 



15.000 



3842.8 



20,000 



4517.1 
4348s 8 
3976.6 



25,000 



35.000 



30,000 



3733.5 



40,000 Service 
ceiling 
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In this case the service ceiling is the optimal 
altitude. But if this altitude can not be reached because 
of poor engine performance in a particular airplane, then 
20,000 fto would be the best altitude. 

An experienced pilot may often fly a near-optimal 
flight profile merely by an educated guess or by carefully 
managing the fuel flow. But he can never do better than 
the computer prediction, and will probably do worse most 
of the time since his wind information is sketchy if he 
has any at all. 

The professional doubter may complain that even mete 
ology doesn't know the accurate wind information at altitu 
so vjhy even bother with a computer solution? We can only 
reply that some good estimate is better than none at all, 
and that modern meteorological equipment can measure and 
predict the winds quite accurately. 

The incorporation of a computer system to predict 
optimal flight profiles may at first meet with some inert! 
from fleet pilots, especially the more experienced ones. 
The "seat of the pants" pilot may reject the computer 
decision as hocus-pocus or just plain Incorrect. But even 
the most experienced aviator may encounter vertigo some 
dark, rainy night and he may discover that flying the 
airplane is about the only job he can handle. Figuring ou 
an optimal flight profile to the beach may take a back 
seat to surviva.1. it is in this situation that informatio] 
from departure control would be most welcome. 
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The cost involved in establishing and operating a 
system such as described above may appear to be too large 
when compared to the small savings in relatively inexpensive 
jet fuelo But the possible loss of a Phantom II and/or 
a pilot may shift the balance in the system's favor. 

One possible reaction against the use of computer 
flight planning for cross-country flights is that often the 
air traffic control center will not allow a climb to the 
desired altitude because they must wait until they can safely 
fit the aircraft into the West-East (or vice-versa) traffic 
flow, which may be considerable in this jet age. This will 
result in remaining at a non-optlmal altitude for many 
minutes. But this drawback exists for non-computer flights 
as well. And it is still better to climb from a non-optimal 
altitude to an optimal one when cleared to climb, than to 
go from one non-optlmal altitude to another. 

Another great advantage in long range computer flight 
planning occurs vjhen several alternate routes are available. 
Then the program can be written to choose the optimal 
altitude as well as the best route of flight. The many 
possible routes and altitudes would require too much compu- 
tation for a pilot during his flight briefing, but a high 
speed computer can solve such a problem In minutes. Thus 
the pilot Can spend his pre-flight time briefing the mission 
and tactics, and let the machine do the arltivnetic. 
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APPENDIX A 



Aircraft performance Charts for the Navy a4c Aircraft 
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MODEL; A-4A, A-4B, A-4C 
ENGINE; J65-W-16A 



DATA AS OF: 1 DECEMBER 1965 

DATA BASIS; FLIGHT TEST (NAVYI 




Figure 6. Climb Fuel 
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PRESSUlfi ALTITUDE - ICOO FEET 



MODEL: A-4A, A-4B, A-4C 
ENGINE: J65-W-16A 



DATA AS OF: 1 DECEMBER 1965 

DATA BASIS: FLIGHT TEST (NAVT) 




CALIBRATED AIRSPEED - KNOTS 



Figure 7 . fillmb Speed Schedule 
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TIME - MINUTES 





Figure 8. Climb Time 
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TRUE MACH NUMBER AVERAGE GROSS WEIGHT- 1000 POUNDS 



MAXIMUM RANGE CRUISE 

TIME AND SPEED 



MODEL: A-4A, A-4B, A-4C 
ENGINE: J65-W-16A 



DATA AS OF: 1 DECEMBER 1965 

DATA BASIS; FLIGHT TEST (NAVY) 




OUTSIDE 

AIR TEMPERATURE- MAXIMUM RANGE TRUE AIRSPEED - KNOTS 
° CENTIGRADE 



figure 9. Maximum Range Cruise - Time and Speed 
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NAUTICAl MILES PER POUND Of FUEL AVERAGE GROSS WEIGHT - 1000 POUNDS 



I 



MAXIMUM RANGE CRUISE — FUEL 



MvOOEL: A-4A, A-4B, A-4C 
ENGINE: J65-W06A 




DATA AS OF: 1 DECEMBER 1065 

DATA BASIS: FLIGHT TEST (NAVY) 




FUEL FLOW - 1000 POUNDS/HOUR 



Figure 10. Maximum Range Cruise - Fuel 
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DISTANCE - NAUTICAL MILES GROSS WEIGHT - 1000 POUNDS 



MODEL: A-4A, A-4B, A>4C 
ENGINE: J65-W-16A 



DATA AS OF: 1 DECEMBER 1965 

DATA BASIS: FLIGHT TEST (NAVY) 



22 



20 



18 



16 



14 



12 



10 







GROSS WEIGHT— 1000 POUl 


SIDS 




DRAG INDEX 


10 


12 


14 


16 


18 


20 


22 


0 


180 


200 


215 


230 


240 


255 


270 


100 


165 


180 


195 


210 


220 


235 


245 


200 


155 


170 


180 


195 


205 


220 


230 



Figure 11 



Descent Distance 

h9 



TIME -MINUTES GROSS WEIGHT - 1000 POUNDS 



MODEL: A-4A, A-4B, A-4C 
ENGINE: J65-W-16A 



DATA AS OF: 1 DECEMBER 1965 

DATA BASIS; FLIGHT TEST (NAVY) 



22 



20 



18 



16 



14 



12 



10 





Figure 12, Descent Time 
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FUEL -POUNDS GROSS WEIGHT _ 1000 POUNDS 



V 



MODEL: A-4A, A-4B, A-4C DATA AS OF: 1 DECEMBER 1965 

ENGINE: J65-W-16A DATA BASIS: FLIGHT TEST (NAVY) 




Figure 13. Descent Fuel 
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APPENDIX B 



Program BINGO and Sample Output 
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r 



BTNGO 



WNnfPS(4''), WNf)'*, ort( -V" ) , riTMn(^.0> 

rjiMFMCfDN HOP! I , TESC(43), Qo AL T ( 4 ": ) t nij T ( 4''' J 

COMMON WNnr.(?S,WNnSPO,CRS 

FTRST data CM?0 : RFAOING AND OISTANCF TT F I El 0 
F'IFL on ROArt, /^trcRAFT type 



P = AD( = ,9'^) C°S ,DI ST ,niFL, ACTVPF 
WPTTC((S,Q0) CPF,OIsf ,FHF| ,ACTY0F 
o/' r0PMAT(4Pi I 
no 9? 1 = 1 , 40, s 

l = M-l 
K = T + '^ 

L = T+R 
M= T +4 

OF An( F , 91 I WNOroc ( 1) t WNO$on( T ) ,WNnr.PS ( J ) , WNDSPO (J I , 
•<^UNnCPS(K) ,WN95Pn(K» ,WNnCPS(H ,WN9S0n(l ) , 
*WNnC‘^S(W),WNDSPD(M) 

<31 format ( 10F<=: , 1 I 
oo roNTTNljP 



PACH FAOn CONTAINS 5 ' TO ppft HF N I NO TNFH 

I P ( AfTYOF. eg. 1 . ) GO TO i^r 
GO TP 41 
4'- RASir = no^''. 

LEVFi. AT irC"- ft 



4 1 



4 0 



OOl 

73? 



73? 



7?4 



A L T = ’ o -> o 

WT=FtlFL*RASlC 
APIJFI =FU'^I 
API ST=ni FT 
aoi ST=ni ST 
T =’ 

R'lRNl I ) =o 
RFO 1 FT = ROT PT 
r^'v^T=WT 



JT=1 
awt=wt 
CONTI MUF 

START= AVHO WT( WT , I , OIF T , Ml 

WT=STAPT 
CAS=HRSPn(WT,T ) 

TAS = CTAS(rAS, T » 

GSPD=WIN0( I ,TAS » 

F1JCL=HRFUEI ( WT , II 
WOIST=REniFT 
TI MF=wni ST/GSPn 
WOI ST = TI MP4-TAS 

I = (WNnCRS(n GE.lR'^.l 00 TO 731 
CMFC<=WNnCPS( 1 » + lR?. 

GO TO 732 

CHECK'^WNOr P S( n-lR"^. 

CONTINMF 

ni FF = CRS-CHPCK 

CORR=ARS(OT FF» 

TF(CPPR,LT 90 .) GO TP ?33 
1C(Corr,GT. ?70. ) GO TO 733 
GO TO 734 
TFST=WDI ST-ni ST 
tf<;t = ars( tf<;t) 

I=( TFST.LF. ? ) GO TO 797 
wt=rwt 

ni ST=ni ST-?. 

GO TO 42 
TFST=W0I ST-OI ST 
TEST=ARS(TF$T) 

I=( TFST.LF 2 ) on TO 797 
WT=RWT 
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ni ST=ni + 
r,n rn 4? 

■^07 cnNTINIJF 

IP( JT .r,T. ’ ) r,n jn 6 3 
H0RT7n )=r.AS 

CHN-timije 

WT=BWT 

SHDI <:T = FnFl *5^ * . 

RF0I<:T=OEni ST-<;Hr>I ST 
APG=AFS ( PFHI ST) 

II =2 

TF(RFni ST> I. T, SHDIST) GO TO 61 
RIIPN( I ) =RHPN( I )+'^ 

RWT^RWT-E')'^, 

WT=RWT 
ni ST=PFniST 
wni ST=PFni«-T 
GO TO 4? 

6 1 BMRN( I ) =RljRN( I ) . ♦ ( ARG /F'FL ) 

Cf IMR( T )=^ 

nesc( I 

WRITF(6,17) 

17 FDPMat( 1H! ,35X, • BURN GLI^R Hnny/ •, 

*RX»* HFSr AIT M 

WP I TP (6 ) 

7C7 POPMA7( 35X,‘ •, 

*^x , • *T 

C 

WP ITF (6, 7? ) B'JRNI 1 ) ,Pl IMP ( 1 ) ,HPR I 7 ( I ) , n ESf ( 1 ) , I 
17 PPPMAT( ////,7RX,4F1.0 4,17) 

r 

r 

r START CLI’^RING IN 10'' i FT tnGRFMFNTS 

r 

>;^t = API)FL + RA SIC 

N=1 

MN = 1 

MM=ir 

on 10'''' 1=2 t'A" 

R'IRNI I )=0. 

JI=1 

CIOIST=0. 

MT = I /MM 

TP(mi.F 0.1 ) GO TO R6 
GO TO 57 
Pf. mm = MM + 5 

r 

C COMPIITF THP CLIMR FUEL AND DISTANCE 

r 

DO 82’-^ J = B,I,5 
ALT = I:kl007 
OPALT ( I ) =AI T 
CAS=C.ISPD( J ) 

TAS=CTAS(CASt J) 

TIMF=CLTIMF(WT, J) 

TIME=TIMF/( J/S) 

TI mE=TIMf/40. 

VFRT=5000. /TIME 
V=RT=VPRT/60RO. 

TAS = S0R7(TAS**2-VPRT**2 ) 

GSPD=AVWIND( J,MN) 
mnj = mn+5 
GSPD=GSPn*7 AS 
F'IEL=CLFUFL(WT,I ) 

RlIRNI T ) =FUEL 
rLDIST=CLDl ST+GSPD*TI me 
R 2r CONTINUE 
R7 CONTINUE 
OIST=BniST 
WT=AFUFL+RASIC 
ALT = I*100''. 

OPALTd )=ALT 
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r'is=rL'^pn( I ) 
rLi«R( n =ffts 

0ELTA = ALT-1 n'v't, 

TAS = r.TAS(CAS, I ) 

TIMExTITT*^*-! WT,I ) 
tt MF = TI 

VPRT=nFLTA/TTMF 
V^^RT = VFOT/<^^^'^, 

T AS = SQR T( TA S^*-'>-UFPT- -7 ) 

r,c;pn=&VWIMn » T , I ) 
r,cpo=FFpn + TAS 
pi)PL=n.FnFL ( wT , n 
BURN( I ) =F'IFL 
CLn!FT=GSPn*TiMF 
C 

C CnvPIlTt^ THF DFSCRNT FIICl ANO OISTAMFR 

r 

TTOI FT=niFT 
AWT = WT-R|)PM ( T ) 

RFFnPF = BI)PM( I ) 

8wt=awT 

WA T T = AWT 
^ 4 ''l = l 
NK=5 

0^ SS*^ J = 1 ,N»^ 

WT=AWT 

r.AS xnSSPOfWAIT) 

DFFC( I) =r AF 
OSniFT = Df^IFT( I ) 

T \S = CTAS(CAS» I ) 
timf=pfdIFT/taf 
MN = l 

^,FPn = AVVITMn ( I ,M\i) 
r,spn=GSPn+TAF 
nSf^T FTsiGFPn + T I MF 
GUT( I ) xDSni ST 
r 

r COMPUTE THF horizontal FOCL AND OI<;T^MCF 

r 

ni ST = TODrST-(Cl. DI ST + f^SniCT ) 

AniST=ni FT 
R'=DIFT = OIFT 

T F ( DT FT, GF* 0. > C,n TO O'; 
oz wRTTF(/S,04> I 

04 PORMATf /,! 2X, • FOP TM!C niFT/FlJEL FOMniNATlOM, *, 
*• stay RELOW ',1'^,' TOnUSANO RFFT •) 

GO TP \ Qni^ 

CONTINUE 

FTART=AVHRWT( WT, I ,niFT, JT ) 

WT=START 
CAF=HPFPn(WT,I ) 

TAS = r.TAS(CAF,I ) 

GFDD=WTNn( I ,TAS ) 

FIIFL = HRFUFl ( WT, n 
wni ST xADIST 
TI ME=Wni ST/GSPP 
WDIST = TIMF*'TAF 

IF( WMnCRS( I ). GE.IRP, > GO TO 4 7(t 
CHECKxWNnCP F( I > + 18'^. 

GO TO 4P7 

CHECK=WNnCPS( I J-IRP. 

i?7 fontimue 

niFF=CPS-CHECK 

r,ORR=ABS(niFF) 

I^fCnPP.LT. ) GO TO 428 
IFICORR.GT. 27'^. » GO TO 428 
GO TO A2P 

428 TFFT=W0IST-0IST 
TEST=ARS( TFST) 

IF( TFFT.LF. 8. J GO TO 481 

WTxAWT 

NN=NN+1 
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ST=ni ST-P. 

GO TO 05 

4?Q TP5T=W0IST-ni«:T 
TCST^ABSC TP<5T) 

IF( TFST, LF. B. ) GO TO 431 
WT=4WT 

NJN = NN+1 

OIST=niST+P. 

GO TO 95 
40 1 rONTTWIJF 

IF( JT.GT.l ) GO TO 7?7 
HOR I Z ( n =r.AS 
CONTTNUF 
JT=? 

FUFL=WDI ST/FUFL 
WAI T=AWT-FMFl. 

55*= CONTTNUF 

B'fRNC I )=R1ICN( I )4-F(IFL 
WT=WA I T 

F'IEL = nSFt)PL C WT, I ) 

B!1RN( T ) =Rl|RN( I ) + F'IFt 

Tc( BURNJC n LF. RltRM( N) » GO TO 5’ 

GO TCI 5 2 

Cl 1^=1 

*=? rnK)TIMI)F 

WP ITT (6, 71) BUPN( I ) ,CL I MR( I ) ,HOP I7( I ) , OFFT ( I) , I 
71 cnPMAT(35X,4Fl 0.4, !7) 

^ C=^NTINUF 

WRTTF(6,73) 0P\LT(N) 

70 CORMATC ///////// ,41X, ’ OPTIMUM ALTITUDE IS 
*• FFPT •) 

WRTTF(6,04) CLTMP(N) 

7^ FOPMATI /,41 X, • Cl I «B SPFFO !<: ',F4 !,' KMOTS •) 
WPTTF(6,75) HOPTZ(N) 

75 POPMATI / ,41 X, • CRUISE SPFFD IS '.Fft.I,' KMOTS ') 
WRITF(5,R6) OUT(N) 

74 PORMATI /,41 X, ' START DESCENT WHFN OUT •) 

WRITF(ft,77) DFSC(N) 

77 FORMAT! /,41X» ' DFSCFNT SdFFO IS • , F6 I,' KNOTS •) 
WRITE(5,7B) BURN(N) 

78 FORMAT! / ,4! X, * FUEL USED WILL RF • . FQ POUNDS *) 

END 



FUNCTION AVHRWT ! WT , I , OIST ,.) I ) 

C 

C COMPUTE THE AHERAGF HORIZONTAL WEIGHT 

r 

hold=wt 

AL7=i*i ono. 

DO 3033 J.) = l»5 
FUFL=HRFUFL !WT, I ) 

FUPL = DI st/fui=l 
WT=HPI 0-! FUEL/?. ) 

I*"! T.GT, 1 ) GO TO 0530 
WRITE!4,3332) WT , I 

3332 FORMAT! Fio. 3,5X, 17) 

3333 CONTINUE 
avhrwt=wt 

WTsHOLO 

RPTUPN 

PMn 
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r 

r 

r 



FM\r ttom n <;po( i > 

THIS IS Thp ri TMR SPCPD FCHFOIH E 



!F( AIT. lf- ’ 6 ‘^sr. ) r,n rn 

CLSPn = 307. - ( ( alt-265' T. J/17 5';'^ • 

r,G TO 1 1 

IP rLSPn=337, -( AI.T/7650' 1*33 

11 CINTTWHF 
RFTUPN 
cMH 



function n ttmfiwt , n 

time nfedfo to climb Tn AiTiTiine 

ALT = I 

IF ( WT. LF. 1?0PP- I GO rn 76] 

TFI WT. LE. 130OC. I GO to 762 
IF( WT. LF. 14001. I GO “^O 763 
IF ( WT. LE. ■' '=OOP. ) GO TO 764 
!F( WT.LE. 16000. I GO TO 765 
CLTIMF=12, 

GO TO 766 

T61 IF( alt. I'F 20. ) r,n to 76U 
tF(ALT.LE.30 ) GO to ^612 

CLTIMF = 7. !-( {40.-ALT) /lO, )*7 7 

GO TO 766 

7611 CLTIMP = ?.7*( ALT/?-. ) 

GO TO 766 

761? ri.TIMF=4.6-( (30. -ALT! /lO. 1*1 7 
GO TO 766 

T62 IF( Al T. LF. n, ) GO TO 7621 
IF( ALT.LF. 30. ) GO TO 7622 
CLTIMF=8. -I (40 -ALT!/ 10, 1*3.2 
GO TO 766 

762! CLTIMF=( ALT/lo. 1^1 . 2 
GO TO 766 

T6 2? CLTIMF=4. 3-( ( ^0 , - A L T > / 20 . I * 3 6 
GO TO 766 

763 IF( ALT. LF. 1 0. ) GO TO 7631 
IF( Al T. LE. 30 » GO TO 2632 
CLT!MF = 9.-( (40. -Al T )/ *''. K3.P 
GO TO 766 

76?? CLTIMF=5.2-( (?p.-ALT> /?0. >*35 
GO TO 766 

26?1 CLTImf=1.6o{ alt/10. I 
GO TO 766 

766 IF( Al T, LF. 15. ) GO TO 76A1 
IF( AIT.LE.20, ) GO TO 7642 
CLTIMF=10. - ( (40. -ALT) /! 0. 1*4 1 
GO TO 766 

7641 CLTIMF = ?.4* { alt/15. ) 

GO TO 766 

764? CLTIME = 5.o-( (3P.-ALT) /I 5. )*3 6 
GO TO 766 

765 IF( ALT.LE, 15. ) GO TO 765! 

IF( ALT, LF. 25 . ) GO TO 7652 
CI.timf = 1 1 . 4-( (40, -alt ) /5. ) * 3. 2 
GO TO 766 

76F1 CLT IMF = 2. 6*( ALT/15. ) 

GO TO 766 

76 5 2 CLT tme = 8. 2-((35-AlT)/20,)*5.6 

766 CONTINUE 
RFTUPN 
ENO 
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^ 0^0 



PtlNCTION HPFDEI. (WT , I ) 

COMPUTE THE FUEL USEO TN LFVFL FLIGHT 
4LT = T 

TP( WT.l E. 130"''^. ) m TO 
IF(WT.LE.l'fOOO. ) GO TO 651 
!«=( WT.LE. 15^00. ) GO TO 654 
IP( WT.LE. 140*''^, ) GO TO 652 
GO TO 652 

651 IF( ALT.LE. 20. ) GO TO 6511 

HPFUFl =.20-((4'',-ALT>/15. >*.''9 
HRFUFL=HRFUEL--G75 *(wT/14C'T0 ) 

GO TO 657 

6511 HPFUFI =.2'^-n?5,-4LT)/25. )*.i 
HRFUFL=HRFUFL-. P'2 *(WT/14''on ) 

GO TO 657 

6‘=5 Ii=( ALT.LE. 20. ) GO TP 6551 

HRPI|FL = .30-( (40. -ALT) / I 5, )*.ro 
HOFUFL=HRFUEL-, 0 2 5* ( WT/1^^ ->0 ) 

GO TO 657 

6551 HRFUFL = .22-((25,-ALT)/25. )*.1 
HRFUPL = HRFUEL-. 0 7*( WT/1 OOO.G. ) 

GO TO 657 

6<=4 TF( ALT.LE, 20. ) GO TO 6541 

HPPUEL = .7T'-( (4P.-ALT) / I 5. ) * . '' 9 
HPFUEL = HRFUEL-,0 2‘^*(WT/16noO. ) 

GO TO 657 

6‘=4l HRFUFL = .19-( ( ?p. -ALT ) /2 5. )*.l 
HRFUFL=HRFUEL-. 02* ( WT M 5^00, ) 

GO TO 657 

65? TF(ALT.LE. 10, ) GO TO 6521 
TF( A) T.LF. 1 5 ) GO TO 6522 
IP( ALt.lE.^O, ) GO TO 6523 
HPFUFL*. 22 

HRFUFL=HRFUEL-.03 *(WT/16000-) 

GO TO 657 

66?] HRFUFL = .l+r ALT/10. )♦. "4 

HRFUFL = HRFUEI -,0'»*(WT/160 00.) 

GO TO 657 

6 62 2 HRFUFL = . 15-( ( I 6 . - AL T ) /5 . ) + . 01 
HRFUFl =HRFUEL- 02 (WT/ 16^00. ) 

GO TO 657 

6 5? 3 HRFUFl =.21-((36.-ALT)/l5- )*.p6 
HPFUFL=HRFUFL-. 079*(WT/160P0. ) 

6'^7 rONTINUF 

RPTIJRN 

PMO 



FUNCTION nSSPO(WT) 

THIS IS THF oescfnt speeh SCHEOULF 

IP( WT.LE. 1207A. ) 00 TO ?l 

nsspn=i 90. ♦ ( ( wT-1 2 ^ 0 ^ . ) /ao^o )*55 
GO TO 22 
?1 0SSPP=199, 

2? CONTINUE 
RETURN 
ENO 
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F'lN'CTTHM Wf Nn( I ,TAS ) 
niMpMsiPN wNnrRS(4(^) , wMn«: on (4'' \ 
nMMOM wNncR'^ , wNOSPOtCRR 

r 

r. ftmdiitF WINO CnMPONFNTS TO rrujvE°T TAS fMTO r,«:pn 

r, 

TF(WNOCP‘:( T ).GF,lRn, ) r,n TO -^0" 
r.HFCK=wNnrP S ( I )+lRn, 
nn Tn 301 

r.HFCK = WNncPS( 1 >-! Bn. 

} rnNTiM!iF 

PI FF=ro?-rHFCK 

CnRR=ABS ( ni FF ) 

T'^(cnRR,iT. P 3 .) r,n rn 3^0 

IF(cnRR.GT, 27-1. ) r,n TP 74-; 

IF(CnRR, r,E IRA.) GO TO 31-.1 
G3 TO 311 

310J CPR R=FnRp-i 37. 

RAn=(CnRR/l RO. )# 3 , 141 
W1 Nn = TAS-CGS( RAO) A=WNnSDn( I ) 
m TO 111 
311 COR R = rORR-o 3, 

R A0 = ( CTRR/'' RO . ) *3. ^ 41 6 
WT Nn = TAS-R IM( PAD) *WNnRPD( I ) 

GO TD 313 

3-5-' RA0 = (rr)RP/1 RO. )^7, 141 (S 

WINP = TAS+rnR(RAn)*WN|ORPO( I ) 

GO TO 313 
34 ^ rORR = CORP-37'' . 

PAD=(CnRR/lRO. ) 1 ' 3 . 141 A, 

WlNn = TAR+RIN( RAO) ♦WA'DRPOf f ) 

31 ■» COMTTN'jF 
RFTtjRM 
PMO 



C 

C 

r 



FUNCTION CTASICAS, I ) 

THIS rONVFPTS CAR INTO TAS 



41 C 
41 1 
41 ? 

41 3 

414 

41 F 

4] 6 

41 7 
4QQ 



ALT = I 

IF(ALa-.lF. 5 ) 
iFIALT.LE- 10. ) 
IF( ALT. LF, 1 5, ) 
TF( ALT.LF, ? 0 . ) 
IF( AI T. IF. ?S, ) 
TF( Al T. lF. 30. ) 
TF( ALT, LF. 3S. ) 
IF ( ALT, LE. 40. ) 
GO TO 417 



GO TO 410 
GO TO 411 
GO TO 41 2 
GO TO 413 
GO TO 414 
GO TO 41 5 
GO TO 41 4 
GO TO 417 



STGMA=l.+( ALT/5. )*.0373 



GO TO 400 

S1GMA = 1, '^773+1 ( alt-5 )/5.) 



* 



OR 44 



GO TO 490 

S1GNA = 1.143 7+((ALT-K.)/=:. )» 0°63 
GO TO 499 , ^ 

SIGMA = 1.24C 2i-( (ALT-15. ) /5. ) *. 1< 94 
GO TO 499 

RIGMA = 1 , 37+ {( ALT- 30 , ) /5 , )* 133« 



GO TO 499 

ST GMA = i ,49 3 8+((ALT-?5 )/5.)«. 1411 

GO TO 49° 

SI GMA = 1 . 6349+ ( (ALT-3' . ) /5. )*. 1615 

GO TO 490 

SI GMA=1 . 7964+ ( ( ALT-35 . ) /F. ) ♦. 21 91 

rTAS = CAS*S TGMA 

RETURN 

ENO 
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PIJNCTTON CLFU^L ( WT » T ) 

rOMpiiTF the FHFL IlSEr* IN THE Cl IMP 



MT = T 

!•=( WT.LE. 1 2'5T^. ) GO TO 61 
I‘=(WT,lF. 13500, ) GO TO 6? 

IF( WT.LE. 145O0. ) GO TO 63 
IF( WT.LE. 15500. ) r,0 TO 64 

CLFUFL=760. - ( (40 -At t)/40. 

GO TP 60 

61 rLFUEL = 5C''. - ( (40. -AL^)/4r, )**;0O. 
GO TP 69 

6? CLFUFL = 56^, -( ( 40 - ALT » / 40 > *56'" . 

GO TO 69 

6^ CLPUEl =620. -( (40 -ALT>/4^- )*6 2'‘. 
GO TP 69 

64 CL‘=UEl=7'^". -((40 ._^lT)/ 4C. »*7f"r. 
6« TP ( I. I T. ) GO TO 691 
TF( I. I E.2^^ ) GP TP 693 
Tp(I.LT,3'') go TP 69-’ 
CLFUEL=CLFllEL+( ALT/4^ 

GO TP 691 

60? CLPUFL=CLFIJEL+ ( AI T/40 , ) *4 0'). 

GP TP 601 

60’ CLFUFL=CLFUEl.+ ( AlT/4^' )*?0">, 

601 CPNTTMIIE 
30 TURN 
ENP 



PUNCTIDN Oni5T( T ) 

PI STANCE CPVFRFP DURING IDLE DESCENT 
AI.T=I 

IF( ALT. LF. 5. ) GP TO 41 
IP ( AIT. LE- 10. ) GP TP 4? 

TP( ALT. LE. ’ 5. ) GP TP 6? 

IF( ALT, LE. ?0. ) GO TO 44 
IF(ALT.LE. 25, ) GO TP 45 
TF( ALT, LE. 30. ) GO TP 46 
1P( ALT. LE. ’5. ) GO TO 47 
DPI ST =R0,+( (ALT-3 5. )/5. I*’’, 

GP TO 48 

41 DPIST=( ( ALT-1. 1/4. ) *6 
GP TP 48 

4? DPI ST = 6. + ( ( ALT-5, I /5 )*8 
GP TO 48 

43 PPIST=14.+( (ALT-’o. ) /5, )*9 , 

GO TP 48 

44 PDTST=?3.+( (ALT-15, )/5. )*] ). 

GP TP 48 

45 DPTST=33.+( (ALT-20. )/5, )*1 1. 

GO TO 48 

46 DDT ST=44.<-( ( ALT-25. ) /5. )*1 8. 
r.P TO 4 8 

4 7 DPI ST =62. + ( (ALT-30. )'5. )*18. 

48 CONTINUE 
RETURN 
ENP 
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non 



FUNCTION nSFUEL(WT»I ) 

CnwpuTF THE FUEL USED IN THE OFSC.PNT 

41 T=I 

TF ( WT.I E. IPOOvO, ) r,0 Tn 771 
TFCWT.lE, 17000. ) GO TO 112 
IP(WT.lE.140^0. J GO TO 777 
IP( WT.LE. ISOOO, ) GO ’^0 IIW 
IF( WT.LE.I'SO'^0. ) GO TO 775 
05FUFL=?10,*( ALT/40. > 

GO TO 111 

7-^1 0SFUFL = ?B0.*(ALT/4O, ) 

GO TO 11 ~> 

772 0SFl)FL = 260.*( ALT/40. ) 

GO TO 777 

77^ nSFUPL=?40.*( ALT/ 40 . I 
GO TO 777 

•'74 0SFUEL=2?5. ♦( ALT/40. ) 

GO TO 11 ~' 

775 0SFUEL=215. ♦(ALT/40. J 
777 CONTI MUF 

return 

ENO 



FUNCTION HRSPDIWT. I ) 

HORIZONTAL S°FEO SCHEOULF 
ALT = I 

IPfWT.LE. 14000 . ) GO TO 551 
IF( WT.LE. 15 QO 0 . ) GO TO 552 
IF(WT,lE. 16000 . ) GO TO 553 
GO TO 557 

5 fl I>=( ALT.IF.^ 0 . ) GO TO 5511 

HRSpn= 7 ?n.+( (^C.-ALT 1 /lO. »* 10 . 
GO TO 557 

5511 HRSP 0 = 230 . ♦(( 30 .-ALT) / 70 . )* 6 C. 
GO TO 557 

552 IF( ALT.LE. 25. 1 GO TO 5521 
HOSPn= 240 .+( ( 40 .-ALT) / 15 . )^ 15 . 
GO TO 557 

5 c;? 1 hRSPO=255. + ( (25.-ALT) / 25. ) *50. 
GO TO 557 

553 IF( ALT.LE. 25 . ) GO TO 5531 
HRSPD=? 45 .+( ( 40 . -ALT) / 15 . )* 15 . 
GO TO 557 

55^1 HR5PD=26C.+( ( 2 5 . -A LT ) /2 5. ) *50 . 
5^1 CONTINUE 
RETURN 
ENO 
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FUNCTION AVWIND(I,MN) 
niMFN<;iON WNDrQS(40) ,WNDS,Pn(AO > 
rnvMPN WNOCRSf WNOSPOtCRS 

CnwpUTF THF AVFRAGF WIND IN CLIMR ANO DESCFNT 

AVWIND=0. 
no N=MN,T 

TF(WNnCRS(N|,GF.180.) GO TO 40-v 
CHECK=WNnCPS(NI+lB^, 

GO TO 405 

404 CHECK = WNnr.RS(N)-T80, 

405 CONTINUF 

OI«=F = CRS-CHECK 
CORR=ABS(DTFP| 

IF (COPR.LT, 9<^. ) GO TO 407 
!P(CORP.GT, 270. I GO TO 408 
IPICORR.GF, IRO. ) GO TO 409 
GO TO 410 
400 cnRP=CORP-180. 

RA0={C0RR/180. )♦?, 141 6 
ADD=-COS{RAO)+WNDS'’n{ N) 

GO TO 444 

4K CGRR = rORR-Q0. 

9A0=( CHOP/ 180. )*8. 1416 
AOO=-SIN(RADI+WNOSPD{ N) 

GO TO 444 

407 RAO=(CORR/1 80. )*8. 1416 
ADO=COS(RAn)*WNnPPn(N) 
r,n TO 444 

408 C0RR=c0RR-270. 

RAn=ir0RR/180. 1*8. 141 6 
AOD=SIN( RA0)*WN0SO0( N 1 

444 AVWINP=AVWI NO+ADO 
AI = I 

AVW!MD=AVWI ND/AI 

RPTURM 

FNO 
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BURN 



ALT 



1520. 3477 
1436.8254 
1388.4229 
1344. 5466 
1304.7539 
1228. 1814 
1195.5281 
1165. 8674 
1138.9519 
1164.0527 
1259. 1255 
1236,3667 
1216. C037 
1198.0518 
1161.8977 
1048. 8237 
1039.0364 
1030.6504 
1023.5820 
1017.7534 
1356.7603 
1428.3113 
1447. 1089 
1427.0532 
1410.2793 
1208. 3538 
1199.5144 
1194.4165 
1191,4863 
1340.2283 
1349.7646 
1354.3699 
1361,1338 
1369. 8821 
1380.4607 
1388.2974 
1398.0271 
1409. 5244 
1422.6638 
1437.3345 



CLIMB 

0.0 

334.7358 
333.6038 
332.4717 
331. 3396 
330. 2075 
329.0754 
327. 9434 
326.8113 
325.6792 
324. 5471 
323,4150 
322. 2830 
321. 1509 
320.0188 
318. 8867 
317.7546 
316.6226 
315.4905 
314.3584 
313. 2263 
312.0942 
310,9622 
309, noi 
308.6980 
307. 5659 
304,0740 
298.2222 
292.3704 
286. 5183 
280.6665 
274.8147 
268.9629 
263.1111 
257.2590 
251.4074 
245.5556 
239, 7037 
233. 8519 
228.0000 



HORIZ 

287. 9998 

285.9998 

283.9998 

281.9998 

279.9998 

277.9998 
275. 9998 

273.9998 

271.9998 

269. 9998 

267.9998 

265.9998 

263.9998 

261.9998 
260,0000 

257.9998 
256.000C 
254. 0000 

252.0000 

250.0000 

248.0000 

246.0000 

244.0000 

242.0000 

240.0000 

238.0000 

236.0000 
234.000C 

232.0000 

230.0000 

229.0000 
228.000C 

227.0000 

226.0000 

225.0000 

224.0000 

223.0000 
222.000C 

221.0000 
220.0000 



OESC 

0.0 

198.0575 
1«='7,9612 
197. 8650 
197, 7687 
197. 6725 
197,5762 
197. 4800 
197. 3837 
196. 9437 
196.8131 
196. 6825 
196. 5519 
196.4212 
196. 2906 
196. 1600 
196,0294 
195.8987 
195. 7681 
195. 6375 
194,7850 
194. 6200 
194.4550 
194.2900 
194. 1250 
193,9600 
193. 7950 
193. 6300 
193.4650 
192.2687 
192.0694 
191.8700 
191.6706 
191.4712 
ic)l. 2710 
191.0725 
190.8731 
190.6737 
190.4744 
190. 2750 



1 

? 

3 

4 

5 

6 

7 

8 
Q 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
2° 

30 

31 

32 

33 

34 

35 

36 

37 

38 
30 

4C 



OPTIMUM ALTITUDE IS 200C0.0 FEET 
CLIMB SPEED IS 314.4 KNOTS 
CRUISE SPEED IS 250,0 KNOTS 
START DESCENT WHEN 34.4 OUT 
DESCENT SPEED IS 195.6 KNOTS 
FUEL USED WILL BE 1017,8 POUNDS 
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